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Abstract. Fluctuations of the soft X-ray background are investigated using harmonic analysis. A section of the 
ROSAT All-Sky Survey around the north galactic pole is used. The flux distribution is expanded into a set of 
harmonic functions and the power spectrum is determined. Several subsamples of the RASS have been used and 
the spectra for different regions and energies are presented. The effects of the data binning in pixels are assessed 
and taken into account. The spectra of the analyzed samples reflect both small scale effects generated by strong 
discrete sources and the large scale gradients of the XRB distribution. Our results show that the power spectrum 
technique can be effectively used to investigate anisotropy of the XRB at various scales. This statistics will become 
a useful tool in the investigation of various XRB components. 
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1. Introduction 



It is likely that the contribution to the XRB by the 



plasma emission is not limited just to our Galaxy. Cen 



ago by Giaccom ct at. (liJrj^ ) was one ol the first discover- 
' ies in the history of the X-ray astronomy. In the seventies, 
■ the X-ray satellites ARIEL Uand HEAO-1 scanned most 
of the sky. A large degree of isotropy of the X-ray flux 
found in those investigations suggested that the origin of 
the XRB was mainly extragalactic and might therefore be 
of cosmological interest. 

At soft energies (below ~ 4keV) imaging optics of 
X-ray telescopes made possible the analysis of the X- 
ray sky with high sensitivity and good angular resolu- 
tion. Observations made with the EINSTEIN satellite 
have shown that a substantial fraction of the soft X-ray 
background originates in discrete sources, mostly AGNs 
( piacconi ct al. 197S , Tananbaum ct al. 1979| ). A further 
major step in the investigation of the soft XRB was done 



The di Huse A-ray background (AKfcSj detected years fc Qstriker (1999| ) and [Phillips ct al. (2000[ ) argue that a 



using the ROSAT X-ray telescope (Trumper 1990). Deep 
ROSAT pointing observations resolved ~ 80% of the 



XRB at 1 keV ( |Hasinger 1998|, [Krautter 1999| ). Using the 
ROSAT observations Hasinger (1992) has shown that at 



energies 1 keV thermal emission by hot Galactic plasma 
represents a distinct diffuse constituent of the XRB. 



non-negligible fraction of the soft (0.5 — 2keV) XRB flux 
is generated by hot gas in filaments and halos surround- 
ing galaxies and various galaxy structures. Thus, the total 
soft X-ray flux is a mixture of extragalactic radiation and 
emission from the Galaxy. The extragalactic part itself 
is non-homogeneous. Apart from the dominating fraction 
produced by various classes of X-ray sources, some truly 
diffuse component is expected. In order to study individ- 
ual components of the XRB it is necessary to disentangle 
them carefully from the total flux. 

Investigation of the X-ray background could provide 
valuable information on the Large-Scale Structure of the 
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Universe (Barcons 1998). The observed XRB flux repre- 
sents the integrated X-ray emission along the line of sight. 
The bul k of the resolve d soft XRB originates at redshifts 
of 1 — 2 ( Barcons 1998 ). It is likely that in the same epoch 
the largest present-day structures collapsed. Since differ- 
ent cosmological models predict a different description of 
the structure formation, the XRB could provide strong 
constraints on these models. 

One of the basic characteristics of the XRB is the 
angular distribution of the background flux in the sky. 
Measurements of the anisotropy of the XRB at various 
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angular scales arc an effective tool in the study of the 
different XRB components. In the present paper we inves- 
tigate the surface distribution of the XRB using spherical 
harmonic analysis. The harmonic power spectrum is one 
of the basic statistical instruments in various problems of 
extragalactic astrophysics. It has not been widely used in 
the studies of the XRB because it requires high quality ob- 



servational material. Recently Scharf et al. (2000) inves 



tigated the all-sky HEA01-A2 data in the 2-10 keV band 
in order to estimate the XRB dipole. Due to the low res- 
olution of the HEAO 1 data, the angular power spectrum 
of the XRB up to 1=20 was calculated. The statistical 
properties of the XRB distribution have been analyzed in 
the past using various correlation methods. Fluctuations 
of the XRB were s uccessfully measured usi ng the autocor- 
relation function flSoltan et al. 1996, 1999"|) . The relat 1011- 
ships between different classes of extragalactic object and 
the XRB were investigated using cross-correlation tech- 
niques. This method proved to be very effective in the 
analysis of relations between the XRB and normal galax- 
ies (e.g. |Lahav et al. 1993| , Roche et al. 1995 , Soltan et 
al. 199§) clusters of galaxie s ( [Soltan et al. 1996| ), IRAS 



galaxies ( |Miyaji et al. 1994 ), and the cosmic microwave 
background (Kneissl et al. 1997). 

The harmonic power spectrum and the autocorrelation 
function are related (Eqs. (12) and (|l9|)). Both statistics 
provide basically the same information on the anisotropy 
of the observed distribution. However, in the actual appli- 
cations to the specific problems both methods give results 
which are complementary, particularly if the data are sub- 
ject to large uncertainties. 

The autocorrelation function (ACF) of the XRB has 
been measured for separa tions between 0.°4 and ~ 6° 
(fcoltan et al. 1996, 19991 ) using the ROSAT All-Sky 
Survey (RASS). Our present measurements of the power 
spectrum of the XRB are based also on the RASS. A 
comparison of the results obtained by these two methods 
should give a better insight into the nature of the XRB 
fluctuations as well as a better understanding of the errors 
affecting the RASS data. To determine the correlations at 
smaller angular separations the ROSAT pointed observa- 
tions should be used (Soltan et al., in preparation). 

In the paper we collect the formulae of harmonic analy- 
sis and give a comprehensive description of the procedures 
used in our computations. We calculate the power spectra 
of the XRB distribution for several regions of the sky and 
for different energy bands. In the future we plan to use 
these results to investigate the nature of the XRB fluc- 
tuations. The amplitudes and characteristic scales of the 
XRB variations determined by means of harmonic analy- 
sis will be used to evaluate various models of the XRB and 
to specify the contribution of the different components to 
the total background flux. 

The organization of the paper is as follows. In Sect. 2 
we describe briefly the observational material: the RASS. 
Because the RASS has been described in detail in a series 
of papers (see the references in the section below) , we give 
here just the basic information. In Sect. 3 the spherical 



harmonic analysis is described in detail. We list in a con- 
sistent way all the formulae and describe the modification 
of the analysis resulting from the partial sky coverage. A 
method to verify the accuracy of our computations is also 
given. A summary of our results is presented in Sect. 4. We 
conclude our investigations with some remarks concerning 
future work where we plan to discuss the implications of 
the present results for various models of the XRB (Sect. 5). 



2. The data 



The RASS has allowed us for the first time to study an 
unbiased, spatially complete sample of the X-ray back- 
ground. 

During the six month survey ROSAT , using the 
Position Sensitive Proportional Counter (PSPC), scanned 
the sky in great circles including the ecliptic poles, collect- 
ing data within the 0.1 — 2.5 keV band. The plane of the 
great circles precessed by ~ 4' each orbit, following the av- 
erage solar motion. Since the full field of view of the PSPC 
had radius of 57', areas of the sky close to the ecliptic plane 
were visible for at least 30 orbits or ~ 2 days, while the 
ecliptic poles were covered by observations during the en- 
tire survey. The exposure therefore increases significantly 
with ecliptic latitude, from a minimum of ~ 700 s at the 
ecliptic plane to over 50 ks at the poles. This variable ex- 
posure introduces strong variations in the signal-to-noise 
ratio. For a full er des cription of th e RASS see Snowden & 
|Schmitt (1990| ) and |Voges (1992"| ). 

Various effects and components contaminating the cos- 
mic signal contribute to the raw data. Non-cosmic counts 
of the PSPC observations consist mostly of charged parti- 
cles and gamma rays which penetrate the detector, small 
pulse-height events that follow larger events in the de- 
tector and are believed to be caused by negative ion for- 
mation and scattered solar X-ray background produced 
by scattering of solar X-rays in the Earth's atmosphere. 
Other sources of non-cosmic background are auroral X- 
rays and other sporadic events, most of which are proba- 
bly also produced by low-energy charged-particle interac- 
tions with the atmosphere or telescope. They last typically 
for a few minutes to a few hours and are called short-term 
enhancements. The last component of non cosmic contam- 
ination is a long-term enhancement - X-ray background 
of unknown character, that varies on a time scale of few 
days and appears most strongly in the Rl and R2 bands. 
For details see |5nowdeii et al. (1994 ) . 

Great efforts were made to separate all these non- 
cosmic counts from the genuine X-ray background 
( Plucinsky et al. 1993j , Snowden et al. 1994, Snowden et al 



1997, and references therein). The final RASS maps con- 



stitute unique material on the large-scale distribution of 
the XRB. However, we are interested in the subtle effects 
which are highly sensitive to even small imperfections in 
the data. Thus, in spite of the careful elimination of the 
non-cosmic counts, effects of any potential residual con- 
tamination could have significant impact on our analysis. 
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Table 1 


Sample 


I min 


f-max 




n[Sr] 


A 


0° 


360° 


12° - 37° T 


3.18 1 


B 


45° 


280° 


12° - 37° T 


2.04 


C 


45° 


280° 


40° 


1.43 


D 


70° 


250° 


40° 


1.12 



Fig. 1- Section of the RASS centered on the North 



t Minimum galactic latitude varies with the galactic longitude 
(see Fig. 0). 

* After removal of region 'R'. 



of sky marked 'R' in Fig. [l] has not been covered by the 
RASS. A large portion of the north galactic hemisphere 
is dominated by the hot plasma emission in our Galaxy 
(Snowdcn ct al. 1995). Most of this local signal is concen- 
trated at galactic longitudes between 260° and 50°. The 
remaining area is visually free from the foreground galactic 
signal, but some local contribution is also present ( [Soltan 



Galactic Pole in an equal area projection The whole map 



covers 3.19 Sr. The areas defined in Table 1 used in the 
computations are marked. The small region marked 'R' 
has been removed from the data because of insufficient 
exposure time. 



Since most of the non-cosmic background occurs in 
the low and medium ROSAT energy bands, typically be- 
low 1 keV (with the exception of the particle contribution 
which becomes noticeable above 2keV), we used the hard 
part of the RASS data, divided into three energy bands: 
R5, R6 and R7, centered at 0.83, 1.15 and 1.55 keV, re- 
spectively. Due to the moderate spectral resolution of the 
PSPC, all bands, particularly the neighboring ones, over- 
lap considerably. Because of this, the count rates in our 
bands are highly correlated. 

The energy band R6 (0.73-1.56 keV) is regarded as the 
best probe for the diffuse cosmological XRB, because in 
this band the systematic uncertainty induced by the fore- 
ground is minimized. We concentrate our analysis and re- 
sults on R6 band, but investigate systematic effects by 
comparison with the two neighboring bands, particularly 
the R5 band. This band is also low in non-cosmic pho- 
tons and in contamination by charged particles, but, as 
the amplitude of the galactic component relative to the 
extragalactic signal increases drastically towards soft en- 
ergies, the R5 band contains an increased galactic fore- 
ground component compared to the R6 band. Thus the R5 
band should allow us to discriminate between the galactic 
and the extragalactic signals. 

The data from the north galactic hemisphere are used. 
The XRB flux is binned into 12' x 12' pixels composing 
an array of 512 x 512 elements. An equal-area polar pro- 
jection centered on the North Galactic Pole has been ap- 
plied to record the data. Within this area a small patch[] 



1 For the present computations this region is defined as 
(287° < I < 290°, 30° < b < 40°) and (288° < I < 293°, 40° < 
b < 47°). 



et al. 1996, 1999). To investigate effects of the galactic con- 
tamination on our determinations of the power spectrum, 
computations have been performed using various subsets 
of the available data. 

The regions investigated in detail in the present analy- 
sis are defined in Table 1. The total area ('A' in Table 1) 
is highly non- homogeneous because of the galactic compo- 
nent. The spherical harmonic transform for this area has 
been computed mainly as a check of our algorithm (see be- 
low). Area 'B' is the largest fragment of the northern celes- 
tial hemisphere apparently free from strong galactic emis- 
sion. However, the border regions of this area are close to 
prominent emission regions. Area 'B' extends also to low 
galactic latitudes (b < 20° at the "corners" of the RASS 
map) and effects of the galactic obscuration are likely to 
be present. Thus, one could expect some interference of 
galactic and extragalactic signals. Area 'C is restricted to 
the galactic latitudes b > 40° . The smallest area , 'D', ha s 
been extensively analyzed by Boltan et al. (1996, 1999). 



The power spectrum obtained for this region should cor- 
respond to the autocorrelation function obtained in those 
papers. To see the effects on our results of the brightest 
sources, some computations were performed without two 
conspicuous clusters: Coma and Abell 1367. 

3. The power spectrum 

Although the definitions and basic relationships for spher- 
ical harmonics are outlined in a number of astrophysical 
publications, we describe below in some detail the stan- 
dard formulae and all the essential steps of our analysis in 
order to present it in a complete and coherent way. Where 
poss ible, we followed the not ation used by Peebles (1973 ) 
and Hauser fc Peebles (1973 ). 



3.1. Sky coverage 

Spherical harmonic functions constitute a set of orthog- 
onal functions on the sphere. Astrophysical data rarely 
cover a full solid angle of 47r. In the present investigation 
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the useful data are also available for a limited region of 
the sky only. Incomplete sky coverage destroys the orthog- 
onality of the spherical harmonic functions. Consequently, 



the power spectrum coefficients are correlated (Hauser & 



Peebles 1973). There are basically two ways to take this 
effect into account. One can either construct a new set 
of functions which are orthogonal in the defined area of 



the sphere (e.g. Gorski 1994, Tegmark 1996), or fill the 
area where the observational material is missing with syn- 
thetic data (see below). Both methods have their advan- 
tages as well as drawbacks. In the present paper we use 
the latter approach for two reasons. Firstly, we perform 
calculations using data selected in several regions of the 
sky. Different windows imply different sets of orthogonal 
functions, which makes comparison of the results ambigu- 
ous. Secondly, we are interested in XRB fluctuations at 
angular scales significantly smaller than the size of the 
area covered by the data. A section of the harmonic power 
spectrum corresponding to these small scales is adequately 
determined using the available RASS data. 

3.2. Spherical harmonics - the formulae 

The distribution of the XRB flux on the celestial sphere, 
p(Q, </>), is represented as an infinite sum of spherical har- 

*7 



P {e,<t>) 



oo I 



1=0 



(i) 



where the coefficients a™ are the spherical harmonic trans- 
form: 



The RASS covers practically the whole celestial sphere, 
but only relatively small regions of the sky contain the 
data dominated by the extragalactic signal; the remaining 
areas are heavily contaminated by local effects. As a re- 
sult of this, we are forced to carry out the power spectrum 
analysis using the data from a fraction of the sky. To as- 
sess the effects of incomplete sky coverage we introduce, 
following Peebles (1973), the window function, W{9,<j>), 
which describes the area used in the analysis: 



w{e,</>) 



1, inside the accepted region, 
0, outside the accepted region. 



(3) 



Using the window function we define for the whole celestial 
sphere a model distribution of the XRB: 



p(0, </>) = W(0, 0) p(0, 0) + [1 - W(0, 4>)] (p) , 



(4) 



where (p) is the mean background flux averaged over the 
accepted area. Thus, within the window area the model 
distribution is identical to the actual one, while outside 
the window it is constant and equal to the average signal 
computed using the data just from the accepted region. 



The model distribution covers the whole sphere, but the 
cosmic information is limited to the window area. 
The spherical harmonic transform of p, 



p Y t m * dfl . 



(5) 



and p is related to the d's by a formula analogous to 
Eq. (0). Using Eq. (g) we get: 



f Q pY™dn-{ P ) f a Y™*dn, for^O, 



4tt (p), 



for I = 0, 



(6) 



where the integrals extend over the solid angle, fi, of the 
accepted area. The zero degree term, o[j, represents the 
average signal in the data and does not depend on the 
window function. One should note that coefficients af 1 , 
although calculated using the cosmic data from the frac- 
tion of the sky, represent a mathematically well defined 
spherical harmonic transform. In particular, the expected 
values of all the I ^ coefficients a™ are equal to 0. Eq. (§) 
implies that 



P Yrdn) = (p) 



(7) 



where (...) enclosing the integral denote the expectation 
value. Eq. (Q) shows that the 'spherical harmonic trans- 
form' coefficients calculated from a fr action of the s phere 
deviate systematically from zero (see [Peebles (1973) ) for a 
detailed discussion). 

One should expect that at angular scales substantially 
smaller than the characteristic size of the window, the 
fluctuations of the XRB are adequately described by a™. 
In order to find a relationship between a™ and a™, we use 



(2) the angular autocorrelation function (ACF), w(8): 



w(6) 



(p(ui) p(w 2 ))| 

w 



(8) 



where a;, = (f>i) for i = 1, 2 denote points in the celes- 
tial sphere separated by the angle 9 and the average values 
(.. .) are calculat ed for the whole sphere. Using Eq. (B2) 
of Peebles (1973 ) and Eq. (0) one obtains the relationship 
between the average flux correlation and the spherical har- 
monic transform: 



(p(wi)p(w 2 ))| 47r = E E 



1=0 m=-l 



P;(C0S( 
47T 



(9) 



{p)2 + j,a+i Zifl(co6j)) 



where 



21 + 



1 1 



m \2 
I I 



(10) 



is the power spectrum for the distribution on the sphere, 
and Pi is a Legendre function of degree I. In an analogous 
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way we define the power spectrum of the model distribu- 
tion p: 



21 



1 1 

— y 
+i ^ 



rn— — l 

Eqs. |, | and M give 



w(9) 



(p) 



-y 

,\2 



i=i 



21 + 1 

47T 



Pi (cos ( 



(11) 



(12) 



and an analogous relationship between w and Zi . If we now 
limit our analysis to separation angles small in comparison 
to the size of the area of solid angle 57, we may express 
the average product of the model flux (p(ui) p(u>2)} \4n as 
the sum of two components (using Eq. (^)): 



4tt- Q 

47T 



(P) 2 ,(13) 



where on the left hand side the averaging is over the whole 
sphere, while on the right hand side it is done over a solid 
angle of fl. One can transform this equation into: 



w{6) 



47T 



(14) 



where wu is the autocorrelation function calculated using 
the data within the solid angle O. As long as the cosmic 
signal contained within the window is a fair sample of the 
entire XRB, wq(9) = w(9). Thus, scaling of the ACFs 
(Eq. (HJ)) together with Eqs. © and @ gives 



4tt 1 



Z, = — 



n 21 + 1 



E 

m— — l 



I5H 2 . 



(15) 



in agreement with Eqs. (53) and (54) by Peebles (1973| ). 
Eq. ([[5]) shows that the high I part of the power spectrum 
(for the whole sphere) is directly proportional to the power 
spectrum based on the data from the fraction of the sphere 
if these data are representative for the whole sphere (fair 
sample assumption). 

3.3. Spherical harmonics - computations 

The large survey area combined with a relatively high an- 
gular resolution of the RASS allows us to calculate the 
power spectrum over a wide range of the spherical har- 
monic degree I. Fluctuations of the XRB on the scale 6 will 
produce a signal in the power spectrum at I ~ tt/9. Since 
the pixel size of the RASS is 12', our algorithm should 
be effective up to I ~ 900. Fortunately this requirement 
could be satisfied with the present day computers without 
major problems. 

The spherical harmonic function is defined in a stan- 
dard way by: 



Y l m {9, 4>) = cf P/"(cos 9) e" n * , 



(16) 



21 + 1 (l-m)\ 



47r (l + m)\ 



1/2 



2 ••: ^^y ^/C^^^ 



o 300 - ... I 



ias = laooj 



300 - 
200 - 

100 -V 



Galactic latitude [deg] 



Fig. 2. Count rate distribution in the R6 energy band 
of the RASS in 12' pixels along galactic longitude I = 
0. The data in this region are heavily contaminated by 
galactic emission which causes large-scale variations in the 
count rate. The curves show the synthesized distribution 
for Z max = 450, 900, and 1800. 



300 
200 




Galaciic latitude [deg] 

Fig. 3. Same as Fig. ||for galactic longitude I = 90°. Most 
of the signal is of extragalactic origin. Because of the negli- 
gible galactic contribution, the average count rate is lower 
than in Fig. ||. Note the good fit for I rmma x — 1800 (bot- 
tom panel). 

P^{ x) = {-i r ^-^p r{x) . 

(I + my. 

where P™ are the associated Legendre polynomials. For 
positive m these arc defined by the formula: 



P ^{x)^{-ir(i~x^—p l {x). 



(17) 



Using the above relations it can be seen that 

>) = {-l) m Yr*{9,4>)- (18) 



I I 
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In the actual computations it is convenient to use recur- 
sive formulae rather than Eqs. (|l6|) and (|T7j) d irectly. We 
use the relations given by Brauthwaite (1972) as quoted 
by Hauser & Peebles (1973). To verify the accuracy of our 



computations and demonstrate the effects of the maxi- 
mum harmonic order on the resolution and quality of the 
fits, we synthesized the model distribution using Eq. (Q) 
and the spherical harmonic transform coefficients a™. 
Because of the finite number of spherical harmonics used 
in the calculations the synthesized map does not repro- 
duce exactly the original RASS data. 



In Fig. ^| we show a sample distribution of count rates 
in 146 pixels along the line of galactic longitude I = be- 
tween galactic latitudes 50° and 80° 0. The three panels 
show the synthesized distributions for three values of the 
maximum spherical harmonic degree used in the calcula- 
tions, Z max = 450, 900, and 1800. The spherical harmonic 
transform a™ was calculated for the total area of roughly 
3.2 sr shown in Fig. |]. In the top panel with l max = 450 
the synthesized distribution follows the large scale fea- 
tures of the real distribution. The minimum angular scale 
'resolved' in this model corresponds roughly to 2 RASS 
pixels. The fit to the real signal looks bad because the 
model locally represents the signal averaged over ~ 4 
pixels and the scan in Fig. [2] has the width of a single 
pixel. The fit is substantially improved in the middle panel 
where l max = 900 and still further in the lower panel for 
^max = 1800. It is evident that to reproduce details of the 
XRB distribution binned into 12' x 12' pixels the harmonic 
decomposition has to reach the degree of l max w 900. The 
model with l max = 1800 reproduces the features of the 
cosmic signal very accurately. We treat this as a verifica- 
tion of our computational algorithm. In Fig. ^| the section 
of the RASS for galactic longitude 90° is shown. Here the 
contamination by the galactic plasma emission is small 
and the signal is dominated by the extragalactic compo- 
nent. The conspicuous large-scale trends are absent in this 
region and the fluctuations are dominated by point-like 
sources. Similar to Fig. || the goodness of the fit improves 
substantially with increasing Z max . In particular, the fits 
for Z max 900 adequately reproduce the peaks associated 
with bright sources. 



3.4. Effects of pixelization of the data 



The power spectrum at high I is affected by the pixeliza- 
tion of the XRB data. The XRB fluctuations at scales 
smaller than the pixel size are smoothed over the pixel 
area modifying the flux correlation at this angular scale 
and removing the variations at smaller scales. We assess 
these effects as follows. 




Fig. 4. A - Power spectrum of all data (Sample A); M - 
model which shows the contribution to the overall spec- 
trum by fluctuations at scales below the pixel size. 

3.4.1. Fluctuations below the pixel size 

The count rate fluctuations recorded in the 12' RASS 
pixels result from the cosmic variations and the photon 
noise. Cosmic variations are produced mainly by discrete 
sources and should be distributed uniformly (apart from 
the source clustering) at high galactic latitudes. On the 
other hand, the amplitude of the fluctuations generated 
by the photon statistics depends on the exposure time 
and varies over the survey area. In principle, this prop- 
erty should help to distinguish between these two sources 
of fluctuations. We plan to analyze this effect in the fu- 
ture paper. Due to the presence of extended X-ray sources 
and source clustering, the cosmic signal binned in pixels is 
correlated, while the fluctuations generated by the photon 
noise are uncorrelated. 

Uncorrelated fluctuations generate the flat power spec- 
trum over all I ^5 900. To assess the contribution to the 
power spectrum by fluctuations at scales below the pixel 
size, we calculate the model spectrum assuming that the 
count rates in pixels are uncorrelated. We use a formula 
analogous to Eq. (22) of |Pccbles (1973| ): 



m |2\ 



w(6) Pi (cos 9) d cost 



(19) 



which is obtained from Eq. (0) above and Eq. (B2) of 



Peebles (1973). It is straightforward to calculate the ex- 



act shape of the autocorrelation function assuming a 'flat' 
X-ray signal within each pixel and no correlation between 
pixels. Such a synthetic autocorrelation function depends 
only on the second moment of the count rates in the pixels 
and on the pixel size (and shape) ; this model autocorrela- 
tion function is equal to zero for separations which are not 
contained within a single pixel. The solid curve in Fig. ^ 

2 Strictly speaking, one side of the pixel row is located ex- 
actly at £ — 0, while the pixel centers are shifted by 6'. 
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Fig. 5. The raw power spectrum of all data. Sample A: 
dotted curve; the corrected spectrum: solid curve; the cor- 
rection for the finite pixel size: dashed curve. 



Fig. 6. Power spectra of the XRB in the R6 band. Labels 
'A', 'B', 'C and 'D' denote the RASS regions defined in 
Table 1. 



shows the power spectrum Z\ calculated for the total area 
(Sample A of Table 1) and the dashed curve represents the 
model power spectrum, assuming that the fluctuations in 
12' pixels are uncorrelated. The observed power spectrum 
is fitted well by the model for I <; 700 and has a system- 
atically larger amplitude for smaller I. This demonstrates 
that the cosmic signal is strongly correlated over a wide 
range of separations. Additionally, good agreement in the 
high-Z range, where potentially numerical problems could 
develop, indicates that our algorithm is accurate over the 
entire range of I. 

3.4.2. Corrections for moderate harmonic scales 

The dashed curve in Fig. |J shows that binning the data 
into 12' pixels reduces the amplitude of the power spec- 
trum also at degrees below I = 900. To account for this 
effect we have calculated the distortion of the power spec- 
trum generated by the binning. The model power spec- 
trum Zi = 1 was used as input data. This flat spectrum 
represents a random distribution of point sources. Due to 
the binning of the input flux, the resulting power spec- 
trum has a shape analogous to the model marked 'M' in 
Fig. ^. The ratio of the modified power spectrum to the 
input model is shown in Fig. ^ by a dashed line. To correct 
for this effect one should normalize the raw RASS spectra 
by this curve. In Fig. |^ the dotted curve represents the 
raw spectrum of the XRB in area 'A', and the solid curve 
shows the corrected spectrum. The spurious curvature of 
the spectrum caused by the data binning is removed. This 
procedure is effective up to a harmonic degree of / 900 
and is applied to all subsequent power spectra computa- 
tions. 



4. Results 

In Fig. | the power spectra of the R6 band are shown for 
four regions of the RASS defined in Table 1. 

The spherical transform at low degrees is corrupted by 
the window function. The maximum angular scale effec- 
tively accessible to the analysis is several times smaller 
than the size of the area used in the computations. The 
strong oscillations of Zi for areas 'B' - 'D' visible in Fig. || 
indicate that the power spectrum estimates of any indi- 
vidual degree I are subject to large uncertainties at low I 
values. Assuming that the 'true' shape of the power spec- 
trum is more regular than that in Fig. |6[ one can try to 
fit a smooth function. This would extend our Zi estimates 
down to I w 10. 

Apart from small / variations none of the spectra ex- 
hibits any distinct features in the entire range of degrees. 
The lack of significant spectral structures is also apparent 
in Fig. [| The RASS data have been accumulated in the 
scanning mode of the satellite operation. Thus the raw 
data consist of a large number of ~ 2° wide strips along 
great circles of constant ecliptic longitudes. In a laborious 
procedure (Snowden et al. 1995 and references therein) 
these raw counts have been reduced to create a homoge- 
neous map of the XRB flux. Despite a great effort to re- 
move all the artifacts produced by the data accumulation 
technique, one should be concerned that not all the instru- 
mental effects have been removed from the final maps. Our 
power spectra do not show visible features at degrees cor- 
responding to an angular scale of 2° , which indicates that 
the XRB maps are free from noticeable effects generated 
by the scanning mode. 

The conspicuous difference of the power spectrum at 
low degree I between area 'A' and the remaining areas re- 
sults from strong contamination by galactic emission be- 
tween galactic longitudes 280° and 45°. One should also 
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note that the power spectrum amplitude for area ! D' is 
systematically smaller than for areas 'B' and 'C. This 
difference indicates that the contribution of the hot galac- 
tic plasma to the background is not confined to the region 
between 280° and 45° but extends further away. It is im- 
portant that, apart from differences at small I, all power 
spectra converge at high I, indicating that fluctuations 
at the smallest scales are dominated by the extragalactic 
component. Although the contribution of the hot plasma 
to the spectrum at larger scales is substantially reduced in 
area 'D' in comparison with the remaining areas, it is still 
possible that some signal is generated locally. To investi- 
gate this problem, we now concentrate on the analysis of 
region 'D'. 

4.1. Error estimates 

Our estimates of the power spectrum are subject to two 
sources of errors. Limited counting statistics and instru- 
mental background subtraction produce uncertainties in 
the distributions of the XRB flux in the RASS data. 
Information on the expected rms count rate uncertainty 
for each pixel has been obtained during the RASS pro- 
duction process and is stored in a separate array. We used 
these data to generate 24 simulated RASS maps. For each 
simulation the count rates in all pixels have been random- 
ized using a Gaussian distribution with an average equal 
to the real count rate and a variance taken from the uncer- 
tainty array. The power spectra of all the simulations have 
been computed. The rms scatter between these simulated 
functions is used as an estimate of the la uncertainty of 
the power spectrum of the real data. The corresponding 
error bars are shown in Fig. [?]a for one of the spectra of 
area 'D'. 

A second source of uncertainty is related to the cos- 
mic variance. To account for this effect, we have divided 
area 'D' into 6 sections and computed the power spectra 
for each region separately. The scatter between the spec- 
tra of the 6 regions is used to estimate the uncertainty of 
the spectrum based on the entire area 'D'. The resultant 
error bars are shown in Fig. 0b. Since the average expo- 
sure times in the RASS are rather short, the uncertainties 
in the count rate for individual pixels are relatively high. 
Nevertheless, one should note that the uncertainties in- 
troduced by cosmic variance are larger than the errors 
generated by poor counting statistics. 

4.2. The power spectrum of the R6 band 

The power spectra shown in Fig. ^ are calculated using 
all the XRB flux in the specified area including strong 
sources. Many of these sources produce count rates which 
exceed the average count rate per pixel by a factor of ten 
or more. In Fig. [8| we have plotted the cumulative distri- 
bution of the count ratesfl in the R6 band for area 'D'. 



Fig. 7. Uncertainties of the power spectrum estimates pro- 
duced by photon counting statistics (a) and cosmic vari- 
ance (b). Area 'D' is used. The error bars represent la 
uncertainties. For clarity, the error bars are plotted for 
selected degree I. 




In ROSAT 



'working' units of 10 X PSPC counts s 



Fig. 8. The cumulative distribution of the count rates in 
the pixels in area 'D' in the R6 band. The solid curve gives 
the number of 12' x 12' pixels with a count rate above the 
specified value; the dashed curve gives the complementary 
distribution. Linear and log scales are used for low and 
high count rate sections of the distribution. The vertical 
dotted line (which splits the plot) marks the median of 
the distribution. Arrows indicate the cutoffs used in the 
calculations (see text). 



The distribution has a long high count rate tail extending 
to above 28 • 10 3 , while the average of the distribution is 
equal to ~ 80. The highest signal, viz. 28175, in the R6 
band is generat ed by a source associate d with the Seyfcrt 
galaxy IC 3599 (|Komossa fc Bade 1999|) . In 0.1 % of pixels 
the count rate exceeds 904 and in 1.0% it exceeds 256. 
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Fig. 9. The power spectra in area 'D' with different 
thresholds imposed on the count rates: 1: all data (same as 
in Fig. |6|); 2: the pixel with the highest count rate has been 
removed; 3: the pixels with the top 0.1 % count rates have 
been removed; 4: the pixels with the top 1.0% count rates 
have been removed; 5; the pixels with the top 1.0 % count 
rates and with negative count rates have been removed; 
the curves in the lower part represent five simulations (see 
text). 

Such peaks generate a signal in the power spectrum over 
all harmonic degrees. To see what impact these strong 
peaks have on the resulting power spectrum we have cal- 
culated the spectra using the data without the brightest 
sources. This was achieved by modifying the window func- 
tion (Eq. (||)) to exclude pixels with count rates above the 
assumed threshold. The effects of the high count rate tail 
on the power spectra are shown in Fig. [)| 

Since the narrow peaks in the distribution produce a 
signal in the harmonic power spectrum over the entire 
range of degrees I, the removal of the brightest sources 
from the area reduces the amplitude of the power spec- 
trum over a wide range of degrees. As a result of the sub- 
traction of the peaks the shape of the spectrum changes 
significantly. For the whole data (Fig. [|, spectrum No 1) 
it is roughly a power law for I < 900, while spectrum No 
4 which represents the data without the brightest sources 
is well approximated by a broken power law with steep 
slope at degrees below ~ 60 and distinctly flatter slope at 
larger I. 

The effects of uncorrelated fluctuations are shown in 
the lower part of Fig. We have plotted there spectra of 
five simulated distributions. In the simulations the count 
rates in the pixels were completely uncorrelated and the 
distribution of the count rates in the pixels was Gaussian 
with a mean value of 75 and a dispersion of 25. As ex- 
pected, the simulated spectra are flat. The scatter be- 
tween the simulations reflecting the stochastic nature of 
the problem indicates the level of uncertainties involved 
in the spectrum estimates. 



0.1 r 



0.001 r 
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Fig. 10. An example of the decomposition of the observed 
spectrum; upper curve: the power spectrum of the data in 
area 'D' after the removal of 1 % pixels with the highest 
count rates (same as spectrum No 4 in Fig. ||); lower curve: 
the residual spectrum obtained from the upper curve by 
subtraction of the flat spectrum representing randomly 
distributed sources contributing most to the high I part of 
the total spectrum (see text). 

Before we investigate this question in detail, it is worth 
noting that the pixels with the lowest count rates do not 
affect significantly the estimates of the power spectrum. 
The low-count-rate part of the distribution is sensitive to 
errors in the instrumental background subtraction. As a 
result of inaccuracies in this procedure, the XRB signal 
formally assigned to more than 1000 pixels is negative. 
However, negative count rate pixels do not create a long 
tail in the pixel distribution and do not influence the over- 
all power spectrum. This is shown in Fig. ^ by spectra No 
4 and 5. In both cases an upper count rate threshold of 256 
was used, and in spectrum No 5 all pixels with negative 
count rates have also been removed. 

A variety of sources contributing to the XRB pro- 
duce its complex structure. Sources which are spatially 
unrelated (e.g. nearby galaxies and distant AGNs or ex- 
tragalactic sources and the local Galactic emission) cre- 
ate patterns on the celestial sphere with characteristic 
and mutually independent statistical properties. Thus, the 
power spectrum of the total XRB signal is a sum of sev- 
eral distinct components corresponding to separate XRB 
constituents. The present data do not constrain strongly 
amplitudes of individual potential components of the mea- 
sured spectrum. To illustrate this we consider spectrum 
No 4 shown in Fig. ^. The flat high degree part of this 
spectrum indicates that a substantial fraction of the sig- 
nal in this range of I is produced by small scale fluctu- 
ations which are generated mainly by point-like sources. 
The distribution of these sources is not constrained. They 
may be distributed nonuniformly and contribute also to 
the large scale fluctuations described by the low I section 
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Fig. 11. The harmonic power spectra of the RASS energy 
bands R5, R6 and R7; area 'D' is used. One pixel (the same 
in all the bands) with the highest count rate is removed 
from the data. 



of the spectrum or they may be distributed randomly so 
that the entire signal at small I is due to the other sources 
contributing to the XRB. In the framework of the latter 
model, there are two separate populations contributing to 
the XRB. The spectra of both components have been esti- 
mated in the following way. Using the procedure described 
in Sect. 3.4.1 we calculated the model spectrum assuming 
that the count rates in the pixels are completely uncor- 
rected. This spectrum was then subtracted from the No 
4 spectrum; the residual spectrum is shown in Fig. |l(| 
The decomposition of the observed spectrum represents a 
model in which strong sources produce local fluctuations, 
while most of the signal at large angular scales is gener- 
ated by a smooth distribution. Its contribution to the local 
fluctuations is negligible. 



Fig. 12. Same as Fig. |TT|; the pixels with the top 0.1% 
count rates are removed. 



4.3. Power spectra vs. energy 

Because the RASS data useful for the extragalactic anal- 
ysis cover a narrow range of energies and the ROSAT en- 
ergy bands overlap substantially, the power spectra of R5, 
R6 and R7 share many common features. Nevertheless 
they exhibit some differences which can be ascribed to a 
slightly different content of the XRB at these energies. In 
Figs. |ll], |l2] and [[3] the power spectra of area 'D' of the 
RASS are shown in three energy bands. The spectra plot- 
ted in Fig. [ll] are calculated after the removal of the pixel 
with the highest count rate in the field (produced by a 
Seyfert galaxy, IC3599); in Figs, [l^ and |l3| the pixels of 
the top 0.1 % and 1.0 % count rates were removed. 

In each subsample the spectra of all three bands, R5, 
R6 and R7 exhibit a similar slope at degrees above I ~ 
200 — 300. Also the normalization of all spectra at high I 
end is consistent with the conjecture that the amplitudes 
of the fluctuations at small angular scales are proportional 
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Fig. 13. Same as Fig. [Lj 
count rates are removed. 



the pixels with the top 1.0% 



to the average signal in each energy band. This implies 
that small-scale fluctuations have the same "colors" as 
the average XRB. 

The similarity of the R6 and R7 spectra continues to 
small degrees I and both these spectra differ significantly 
from the spectrum of the softest band R5. The latter spec- 
trum exhibits a sharp increase below / ~ 20 in Figs, [ll] and 
I ~ 100 in Figs. [l2| and |l^. This is a result of the contam- 
ination of the R5 band by soft photons, which evidently 
exhibit large-scale non-uniformities. The spectra of sub- 
samples devoid of the strongest sources reveal a change 
of slope at quite small angular scales of a few degrees. 
Evidently, the excess fluctuations of the soft XRB over 
the harder component represent a separate XRB compo- 
nent whose origin deserves careful examination. 
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Fig. 14. The autocorrelations functions of area 'D'. Labels 
'2' and '4' denote the ACFs computed from the power 
spectra in Fig. || marked '2' and '4', respectively. The ACF 
calculated directly from the data in subsample No '2' is 
shown by open circles. 



4.4. The autocorrelation function 

The power spectrum may be used to calculate the auto- 
correlation function (ACF). In the previous sections the 
relationship between those two statistics has been used to 
elucidate the effects of the data binning. Here we employ 
this relationship (Eq. (p^)) to calculate the ACF for the 
whole separation range and to compare it with the ACF 
determinations given by {Soltan et al. (1996, 1999 ). The 
relationship between ACF and harmonic power spectrum 
is discussed in the Appendix. 

The ACF computed using two power spectra, No 2 
and 4 in Fig. ^| is shown in Fig. |l4|. For comparison, open 
circles show the ACF computed directly from the RASS 
data using the same subsample as for the curve labelled 
'2'. It is clearly visible that the removal of the strongest 
sources from the area reduces significantly the ACF at 
separations below ~ 1°. This implies that these sources 
are not distributed randomly. The residual ACF signal 
shows that the XRB devoid of the brightest sources still 
exhibits fluctuations. But the substantial decrease of the 
ACF amplitude after the removal of 1 % of pixels indicates 
that the origin of the XRB fluctuations requires further 
studies. 



5. Prospects for the future 

Deciphering the XRB fluctuations plays a crucial role 
in the investigation of various components of the back- 
ground. In the present work we have calculated the har- 
monic power spectrum of several subsamples of the RASS. 
The data have been selected according to various criteria 
such as the energy band, amount of the galactic contami- 
nation and presence of strong sources. The corresponding 



power spectra reveal the effects of the different selection 
criteria. This gives good prospects to use the spectra to 
investigate models of the XRB. 

Such models should specify the relative contributions 
of various classes of discrete sources contributing to the 
XRB as well as the amount of the thermal emission gen- 
erated by hot gas which according to Cen & Ostriker 
(1999| ) is accumulating in the gravitational potential wells 
of galaxies and galaxy agglomerations. Since the tempera- 
ture of the gas is relatively low one expects that the energy 
spectra of discrete sources and of thermal emission are 
distinctly different. Also a smooth gas distribution would 
not produce sharp peaks in the XRB. These factors could 
help to differentiate between the two XRB components 
by means of the power spectrum and ACF analysis. Such 
investigation should also utilize the small scale XRB fluc- 
tuation measurements based on the ROSAT pointings. 

APPENDIX 

The ACF amplitude difference at small separations be- 
tween ACFs in Fig. |lj is in agreement with the expected 
correspondence that the larger amplitude of the spectrum 
at high degree I generates the stronger signal of the ACF 
at small separations. However, the relationship defined by 
Eq. ( |l2| ) does not yield a plain dependence of the ACF 
slope on the slope of the power spectrum Zi. Generally 
the power law shape of the spectrum does not produce 
the power law shape of the ACF. However, for the special 
case of Z\ ~ l^ 1 Eq. fll2|) generates an ACF which is well 
approximated by a power law with a slope of —1. The ex- 
act solution differs from the power law w(9) ~ 1/9 by less 
than 5% for separations 9 < 18? 5. 

One should note that distinctly different harmonic 
power spectra could produce similar ACFs. For instance, 
the ACFs calculated using spectra plotted in Fig. [l(] are 
identical for separations above ~ 0?3 despite large and sys- 
tematic differences in the input spectra. This is because 
the flux correlations between pixels are the same in both 
cases. The lower spectrum has been obtained from the up- 
per one by subtracting the flat spectrum which represents 
random (uncorrelated) fluctuations. Thus, both spectra 
describe distributions with the same correlation proper- 
ties. 
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